High Contrast X-ray Speckle from Atomic-Scale Order in Liquids and Glasses 
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The availability of ultrafast pulses of coherent hard x-rays from the Linac Coherent Light Source opens new 
opportunities for studies of atomic-scale dynamics in amorphous materials. Here we show that single ultrafast 
coherent x-ray pulses can be used to observe the speckle contrast in the high-angle diffraction from liquid Ga 
and glassy Ni2Pd2P and B2O3. We determine the thresholds above which the x-ray pulses disturb the atomic 
arrangements. Furthermore, high contrast speckle is observed in scattering patterns from the glasses integrated 
over many pulses, demonstrating that the source and optics are sufficiently stable for x-ray photon correlation 
spectroscopy studies of dynamics over a wide range of time scales. 

PACS numbers: 41.60.Cr, 61.20.-p, 61.43.Dq, 51.50.th 



Ever since the initial observation of speckle in the diffrac- 
tion of coherent short- wavelength x-rays HI, it has been rec- 
ognized that coherent x-ray techniques have the potential to 
provide powerful new probes of atomic scale structure and 
dynamics in non-crystalline systems, analogous to the tech- 
niques developed to study disorder at much larger length 
scales using coherent visible light |2y. In particular, x-ray 
photon correlation spectroscopy (XPCS) has been developed 
using high-brightness synchrotron x-ray sources to observe 
equilibrium and non-equilibrium dynamics in diverse systems 
IS 01- A limiting factor in the development of XPCS for 
atomic scale studies has been the signal level that can be ob- 
tained. Almost all XPCS studies to date probe relatively large 
length scale (> 10 nm) structures, for example via small-angle 
scattering from colloids and polymers |5j-|9Jl or diffuse scat- 
tering near Bragg peaks of crystals | 1(J- 1211 because of their 
larger scattering cross sections compared with atomic-scale 
fluctuations, and the limited coherent x-ray power available. 
The first XPCS observation of truly atomic-scale dynamics 
was obtained only recently in a study of diffusion in Cu- 
Au on relatively slow timescales (> 10 s). 

One of the scientific drivers for the new generation of in- 
tense, coherent hard x-ray free electron laser sources, such as 
the Linac Coherent Light Source (LCLS) at SLAC National 
Accelerator Laboratory, has been their potential for XPCS 
studies to open a new frontier at the natural time scales of 
even the fastest condensed matter systems, e.g. atomic dif- 
fusion in liquids ffl 151. The combination of femtosec- 
ond pulses and atomic resolution speckle provides unprece- 
dented opportunities to test theories of liquid structure and 
dynamics, such as the intriguing predictions that have recently 
emerged regarding the complex dynamics of both diffusive 
flrill and vibrational [17] modes in liquids and glasses. Ul- 
trafast XPCS provides a time-domain probe complementary 



to inelastic scattering, capable of studies of sub-micron vol- 
umes and non-equilibrium dynamics. To reach femtosecond 
time scales, which are much faster than the time resolution 
of imaging detectors, a pulse split-and-delay technique has 
been proposed Uil [l5 , 18-2CJ] in which the correlation time 
is determined from the change in contrast of the sum of two 
speckle patterns, as a function of the delay time between the 
patterns. Although LCLS provides a huge leap in available 
coherent flux and accessible time scales, the high number of 
photons delivered in a single pulse poses challenges in avoid- 
ing x-ray effects on the sample during a dynamics measure- 
ment [I4j,LL5[]. A primary issue is whether accurate contrast 
values can be extracted from the relatively weak speckle pat- 
terns expected in optimized XPCS studies. Likewise, the vari- 
ations in pulse energy and position pose potential challenges 
for speckle measurements. Here, we report the first obser- 
vation of high contrast x-ray speckle from the atomic scale 
structure in both liquid and glass samples. These were ob- 
tained using femtosecond hard x-ray pulses from the LCLS. 
We present an effective analysis procedure to extract contrast, 
show the dependence of contrast on experimental conditions, 
report the observed perturbation thresholds of various sam- 
ples, and discuss the design of an optimized atomic resolution 
XPCS experiment. 

We recorded speckle patterns from three samples: Ga 
heated into the liquid phase at 35 C, and two glasses at room 
temperature, Ni2Pd2P and B2O3. The average structure fac- 
tors of these systems have all been previously characterized 
using standard scattering methods l2lU23ll . Experiments were 
carried out with the XPP instrument at LCLS, using a Si (1 1 1) 
monochromator with 1.4 x 10~ 4 bandwidth set to 7.99 keV 
photon energy, and a 70 fs electron pulse width 12411 . The aver- 
age unattenuated x-ray pulse energy was typically ~ 1 x 10 10 
and ~ 4 x 10 10 photons per pulse for single-pulse or multi- 
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pulse images, respectively 12411 . although there was a very 
broad distribution of pulse energies with a maximum typically 
5 times the mean. The scattering geometries were chosen 
to provide coherent illumination conditions that would ful- 
fill the requirements to produce high-contrast speckle patterns 
125, 26| and to resolve this pattern at the detector. Measure- 
ments were made with a beam focused to 1.7 /im FWHM at 
the sample using Be compound refractive lenses. Diffraction 
patterns were recorded using a direct-x-ray-detection CCD 
with 20 [im square pixels at a distance from the sample of 
L = 18 to 142 cm, positioned at an average scattering angle 
corresponding to the peak in the strongest amorphous scatter- 
ing ring. For Ga and Ni2Pd2P, we employed bulk samples in 
reflection geometry at Q ss 2.60 and 2.95 A -1 , respectively; 
for B2O3, we employed a fiber of 18 /im thickness in trans- 
mission geometry at Q s» 1.62 A -1 . Unattenuated single- 
pulse speckle patterns were recorded from all samples. To in- 
vestigate the stability of the experimental setup and perturba- 
tion of the sample by the x-ray beam, we also recorded speckle 
patterns from the glass samples that were sums of multiple at- 
tenuated pulses. 

Figure 1(a) shows a typical region of amorphous scattering 
viewed by the detector at L = 37 cm for liquid Ga. This im- 
age of the average scattering was obtained by summing many 
single-pulse patterns and binning the pixels to smooth out the 
photon statistics. Figure 1(c) gives a small region of one of 
the single-pulse patterns, showing the individual photon hits. 
Since the width of the x-ray pulses is expected to be shorter 
than the time scale of equilibrium diffusive and vibrational 
motions of the atoms in the liquid, such single-pulse diffrac- 
tion patterns should exhibit high speckle contrast if the x-ray 
optics are sufficiently free of aberrations, the detector resolu- 
tion and noise characteristics are adequate, and the x-ray pulse 
does not disturb the atomic positions in the sample on the time 
scale of the measurement. A critical issue is to understand the 
accuracy with which the contrast can be determined from such 
weak speckle patterns. 

Recent studies have obtained the contrast in speckle pat- 
terns from the normalized variance of the signal distribution 
on the detector 1271 12811 . When the photon density is low, 
the variance contains a large contribution from the count- 
ing sta tistics of individual photons in addition to the contrast 
Ji , 18 , 24 ] . Our approach to extracting the contrast is to com- 
pare the observed probability for k photons arriving in a single 
detector pixel to the negative binomial distribution expected 
for a speckle pattern of a given contrast lf29tl 
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where k is the mean photon density (counts per pixel), M 
is the number of modes in the speckle pattern, and T is the 
gamma function. The speckle contrast factor (3 = M has a 
maximum value of unity for a single-mode distribution, and 
it drops to zero as M increases. For M — > 00, Eq. (1) 
approaches the familiar Poisson distribution for uncorrected 
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FIG. 1 : (Color online) (a) Average scattering of liquid Ga integrated 
over 309 pulses from LCLS. Detector at L = 37 cm from sample 
subtends scattering angles including the main amorphous scattering 
ring at Q = 2.60 A -1 . Dashed box shows region of pixels used in 
contrast analysis, (b) Average photon density vs. scattering angle, (c) 
Expansion of small region of a single-pulse scattering pattern show- 
ing signal in analog-to-digital units (ADUs) from individual photons, 
and assignment of photon positions. 

photon positions, P{k) = k k exp(— k)/k\. The contrast fac- 
tor will be less than unity if the conditions are not met for the 
required transverse or longitudinal coherence of the incident 
beam 125112611 . if the detector pixels are too large to fully re- 
solve the speckle, or if the sample is not static on the time 
scale over which the pattern is recorded. It is this latter effect 
that allows the sample dynamics to be extracted from mea- 
surements of speckle contrast. In particular, if we record the 
sum of two speckle patterns of contrast factor /3n from two 
equal-energy x-ray pulses separated by a time r, the contrast 
factor of the sum will drop from /3o to /3q /2 as t is varied from 
much less than to much more than the correlation time of the 
sample dynamics 1 15, 29tl . 



To accurately extract the contrast in these weak speckle 
patterns, we first employ a "droplet algorithm" 124l 1 3 Oil to 
locate the positions of each detected photon by fitting each 
recorded pattern. Typical results for extracted photon posi- 
tions are shown in Fig. 1(c). Experimental probabilities P(k) 
are obtained by choosing a region of pixels that has nearly uni- 
form k and Q, i.e. along the peak of the amorphous scattering 
ring as shown in Fig. 1(a), binning the photon positions into 
the detector pixels, and determining the fraction of pixels that 
have k photons. Figure 2 shows the experimental probabilities 
for k = 1 to 4 for 309 speckle patterns having various mean 
count densities k. We obtained a contrast factor f3 for each 
speckle pattern by analyzing the P(k) as described in the sup- 
plemental material 12411 . The weighted average of these val- 
ues gives (/3) = 0.276 ± 0.004 for liquid Ga and L = 37 cm. 
While the speckle patterns with highest k provide the high- 
est accuracy contrast values, the full distribution of weaker 
and stronger speckle patterns obtained (due to the wide inci- 
dent pulse energy variation characteristic of monochromatic 
hard x-ray experiments at LCLS) verifies the predicted k de- 
pendences given by the negative binomial distribution. These 
measurements on liquid Ga represent the first observation of 
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FIG. 2: (Color online) Observed probability of k — 1 to 4 photons 
within a pixel as a function of mean photon density k, for 309 single- 
pulse speckle patterns from liquid Ga for L — 37 cm. Thick curve 
is negative binomial distribution with (3 — 0.276; thin curves show 
limiting distributions for /? = 1 and ft = 0. 



TABLE I: Experimental parameters, extracted average contrast fac- 
tors (/?), and calculated maximum contrast factor f3 ca i c for speckle 
patterns from Ga liquid, Ni2Pd2P glass, and B2O3 glass. The pa- 
rameter L is the sample-to-detector distance, "pulses" is the number 
of consecutive LCLS pulses integrated in each speckle pattern, n pa tt 
gives the the number of speckle patterns included in the contrast anal- 
ysis, (7o) gives the average number of photons per pulse incident on 
the sample (after any attenuation), and (k) is the average of the mean 
photon densities in this set of patterns. The number of pixels used 
was n pix = 3.6, 7.2, or 13.2 x 10 5 for L = 18, 35 to 37, or 142 
cm, respectively. 
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speckle from the atomic-scale order in a liquid. 

To investigate whether the stability of the source and ex- 
periment will affect the contrast in multi-pulse XPCS exper- 
iments, and to understand the thresholds for irreversible per- 
turbation of the sample structure by the incident x-ray pulses, 
we studied two glasses where the liquid-like atomic arrange- 
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FIG. 3: (Color online) Contrast factor (/3) vs. sample-to-detector 
distance L for the Ni2Pd2P glass. Curve is calculation for a static 
sample; symbols are observed values for single-pulse and 500-pulse 
images with different (Io) values shown in legend. 

ment should be static. For the metallic glass Ni2Pd2P, we 
collected not only single-pulse speckle patterns without at- 
tenuation of the incident beam (like we did with liquid Ga), 
but also sums of 500 consecutive attenuated pulses. The val- 
ues of ((3) extracted for both single-pulse and attenuated 500- 
pulse speckle patterns at three different values of sample-to- 
detector distance L are given in Table I. Scattering from the 
low-atomic-number glass B2O3 was too weak to obtain val- 
ues of (j3) from single pulses under the available experimen- 
tal conditions; Table I gives results for multiple-pulse patterns 
with various (Iq). 

The observed (j3) values can be compared with those cal- 
culated for a static sample using an incident beam with full 
transverse coherence but non-zero photon energy bandwidth, 
and a detector with non-zero pixel size. These resolution ef- 
fects reduce the calculated static contrast factor /3 ca i c from 
unity depending upon the size of the illuminated volume, the 
scattering angle, and the sample-to-detector distance L 112411 . 
Table I gives f3 ca i c values for each of the samples and exper- 
imental conditions. The observed (ft) value for liquid Ga is 
90% of the corresponding j3 ca i c , indicating that the ultrafast 
x-ray pulses have effectively frozen the motion of the atoms, 
giving a coherent diffraction snapshot of their arrangement in 
the liquid. It also indicates that the unattenuated, focused x- 
ray beam does not significantly alter the atomic arrangements 
in liquid Ga during the time of the pulse, even though the en- 
ergy deposited in the illuminated sample volume (~ 30 eV 

per atom per average pulse) is typically enough to vaporize it 

I — I 

on a longer time scale [24]. 

Figure 3 shows the measured contrast factor values for the 
Ni2Pd2P glass compared with /3 CQ ; C as a function of L. The 
overall dependence of the observed (/?} on L agrees with that 
of the calculation. The highly attenuated 500-pulse patterns 
have contrast factors that are 75-80% of /3 CQ ; C , indicating not 
only that there is minimal disturbance of the sample by the x- 
ray pulses, but that the contrast is not greatly lowered by vari- 
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ations in the pulse energy and position of the incident beam 
or other instabilities of the experimental setup. The measured 
(P) values are somewhat smaller for 500-pulse patterns with 
less attenuation or single pulse patterns with no attenuation, 
consistent with some disturbance of the sample structure by 
the incident beam. For the single-pulse data, the energy den- 
sity deposited is higher for Ni2Pd2P (~ 75 eV or ~ 95 eV 
per atom per average pulse) than for Ga owing to its shorter 
absorption length, which agrees with the observed higher dis- 
turbance during single pulses [24]. Based on the maximum 
incident pulse energies in the attenuated 500-pulse sequences, 
the threshold for longer time (inter-pulse) sample disturbance 
effects occurs at about 1 eV per illuminated atom absorbed en- 
ergy. The dataset for B2O3 indicates that this threshold is also 
about 1 eV per atom in this "light" material, corresponding 
to a higher incident fiuence because of its longer absorption 
length. Based on the observed melting threshold for Ni2Pd2P, 
we estimate that the deposited pulse energy is thermalized in 
a volume 50 times larger than the illuminated volume, consis- 
tent with the expected spread of the electron/photon/phonon 
cascade generated by x-ray absorption [24], resulting in a 
maximum temperature rise of 0.02 eV/3fce ~ 80 K. 

These measurements allow an analysis of the feasibility of 
femtosecond XPCS experiments using the split-pulse tech- 
nique 

ElHllll. Thefi 

sure of merit for such measurements 
is the signal-to-noise ratio for f3, which for low k is deter- 
mined by fluctuations in the small number of k — 2 events and 
can be expressed as 024J] /3/cra = Pk[n p i X n pa tt / ^{X + Z?)] 1 ^ 2 - 
This expression agrees well with the observed accuracy of (/?) 
given in Table I, indicating that other potential experimental 
contributions to the uncertainty are negligible. If we use a 
detector capable of recording n p i X = 10 6 pixels at the full 
LCLS repetition rate of 120 Hz, to give n pat t — 4 x 10 4 in 
six minutes per delay time, and extrapolate the measured rela- 
tionships between k and Iq Il24ll . the pulse energies needed 
to give a signal-to-noise ratio of 5 sufficient for XPCS are 
Iq = 6,0.3, or 0.13 x 10 8 photons per pulse, correspond- 
ing to sample temperature rises of 3 to 7 K, for B2O3, Ga, 
or Ni2Pd2P, respectively. These temperature rises are signifi- 
cantly below the 80 K threshold for structural disturbance by 
diffusive atomic rearrangement that we observed for the glass 
samples, indicating that even more sensitive samples and pro- 
cesses can be studied. 

The new analysis technique presented here allows the accu- 
rate determination of speckle contrast using the weak patterns 
obtained under conditions in which the x-ray pulses from a 
free electron laser do not disturb the sample dynamics. We 
observe high contrast factors in x-ray speckle patterns from 
liquids and glasses obtained with ultrafast pulses, demonstrat- 
ing that XPCS can be used to observe their atomic-scale dy- 
namics on times scales down to the femtosecond range. We 
find that the atomic motions associated with x-ray damage oc- 
cur on time scales shorter than the x-ray pulse width used here 
for deposited energy densities greater than about 50 eV per il- 
luminated atom 112411 . This provides guidance for design of 
studies in which single-pulse speckle patterns are analyzed to 



obtain atomic-scale structural information. 
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